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Abstract. In a magnetic field, spin-ladders undergo two zero-temperature phase transitions at the critical 
fields H c i and H C 2- An experimental review of static and dynamical properties of spin-ladders close to 
these critical points is presented. The scaling functions, universal to all quantum critical points in one- 
dimension, are extracted from (a) the thermodynamic quantities (magnetization) and (b) the dynamical 
functions (NMR relaxation). A simple mapping of strongly coupled spin ladders in a magnetic field on the 
exactly solvable XXZ model enables to make detailed fits and gives an overall understanding of a broad 
class of quantum magnets in their gapless phase (between H c \ and H C 2). In this phase, the low temperature 
divergence of the NMR relaxation demonstrates its Luttinger liquid nature as well as the novel quantum 
critical regime at higher temperature. The general behavior close these quantum critical points can be tied 
to known models of quantum magnetism. 

PACS. 75.10Jm Quantized spin models - 75.40.-s Critical-points effects - 76.60.-k Nuclear magnetic 
resonance 



1 Introduction 

It is well known that long-range order is destroyed by 
quantum fluctuations in one-dimensional antiferromagnets. 
If the importance of quantum effects is ubiquitous in one- 
dimension, a wide variety of ground states can neverthe- 
less be found in nature. Some systems have a continuum 
of low energy modes, some have an energy gap above a 
unique ground state, other dimerize. Where do these dif- 
ferences come from? In simple terms, the role of quantum 
effects is simply to "connect" different classical ground 
states (for example the Neel states | |, 4)1\ !>•••) an d 
| |, t) I) T) • • •)) by tunneling processes. Depending of the 
strength of the tunneling matrix elements, which can usu- 
ally be measured by a coupling constant g, the system will 
be more or less localized around a classical ground states. 
As g is varied, the system can delocalize at a critical value 
g c . When the system delocalizes in spin space, the ground 
state becomes a rotationally invariant singlet and in all 
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the cases which will be considered here, an energy gap 
A appears simultaneously in the energy spectrum. Many 
physical aspects determine the strength of quantum fluc- 
tuations. The integer or half-integer nature of the spin con- 
sidered modify drastically selections rules for quantum- 
processes [[!]. This is why integer-spin chains for which g 
exceed g c have an energy gap, while half-integer spin- 
chains remain gapless {g < g c ). Other physical parameters 
(exchange constants, applied magnetic fields) also enter in 
the precise determination of the coupling strength g. Sys- 
tems for which the coupling constant can be continuously 
varied by an experimentally controllable parameter, such 
as a magnetic field, are rare. In this paper we review a few 
1-D antiferromagnets where such zero temperature criti- 
cal points have been observed, with a particular emphasis 
on Cu 2 (C5Hi 2 N 2 ) 2 Cl4 (also known as CuHpCl)|,|,|, a 
spin-ladder compound, where a complete set of experi- 
ments exist. 

At a quantum critical point the system 

switches from one ground state into another. Specifically 
when g < g c , antiferromagnetic correlation functions de- 
cay as power laws and the system is nearly ordered. When 
g is increased above g c , a gap opens up and the range of 
antiferromagnetic correlation become finite, of the order of 
£ g ~ a/\g — g c \ (a is the lattice constant). In the vicinity of 
g c , one has to go to relatively long lengthscale, exceeding 



Fig. 1. Phase diagram of spin-ladders in a magnetic field. The 
magnetic field (vertical axis), can be thought as a tuning pa- 
rameter for the quantum coupling constant g. Below H c i g c , 
the coupling constant g exceeds g c (H < H c i <4> g > g c ). This 
phase has a singlet ground state and an energy gap. Above 
H c i, the coupling constant g drops below g c : the phase is mag- 
netic and gapless and belongs to the same universality class 
as the Heisenberg XXZ model. The dotted lines represent the 
crossover lines where £ s ~ £t separating two quantum regimes 
just described (where £ 3 is the relevant quantum correlation 
length) from the quantum critical regime dominated by £t- For 
the system discussed in this paper the quantum critical behav- 
ior at the upper critical field H C 2 is similar and will be discussed 
in the next sections. The shaded area between critical fields is 
a 3D ordered phase stabilized by transverse interactions (see 
Sec. 5). 



£ g to be able to tell in which phase the system is. In other 
words, the nature of the ground state is manifest only 
at long lengthscale. At finite temperature spin-flip pro- 
cesses become possible. They cut spin-correlations off at 
a lengthscale £t which can be estimated in the quantum- 
disordered phase (g > g c ) as the mean distance between 
excitations. Since their energies are higher than the energy 
gap A ~ 1/ Is — 3c | above the ground state, their density 
is activated. Hence when k-sT < A, the mean distance be- 
tween excitations £t greatly exceed £ s and thermal fluc- 
tuations are not really relevant. On the other hand, when 
k^T > A, their density is governed by the relative value 
of k^T — A compared to the bandwidth of the triplet ex- 
citations. When A is small compared to k-gT and this 
bandwidth, £ g exceeds very rapidly £y. In this case, the 
density of excitations are determined by k^T alone which 
become the only relevant energy scale. In this limit, dy- 
namical properties are similar to those of a simple para- 
magnet (1/Ti — > est) while thermodynamics quantities 
remain nontrivial. This regime is (improperly) named the 
quantum critical regime, because most properties are de- 
termined by the single lengthscale £t as in ordinary phase 
transitions. 



To summarize, there are two-relevant lengthscales at 
a quantum critical point, the quantum correlation length 
£, g and the thermal length £j\ Depending on their rela- 
tive values, different regimes exists. They are represented 
graphically in Fig. [j], on the H-T phase diagram appropri- 
ate to spin-ladders. The regions dominated by quantum 
effects are (a) the gapped spin-liquid phase below line 
A, (b) the XXZ or Luttinger liquid phase to the left of 
line C and (c) the gapped polarized phase above line B. 
The quantum critical region is found to the right of these 
crossover line and extends down to T — at the critical 
fields H c i and H c2 . Each regime will clearly be identified 
using thermodynamic and 1 jT\ NMR relaxation measure- 
ments which allow to place precisely the crossover lines on 
this phase diagram. 

The differents sections are organized as follows: a brief 
description of several families of gapped antiferromagnets 
having an H-T phase diagram similar to the one shown in 
Fig. 1 can be found in Sec. 2. A description of the struc- 
ture and the interactions relevant to CuHpCl, the ID lad- 
der system chosen for our case study, is given in Sec. 3. 
A mapping of strongly coupled ladder in a magnetic field 
onto the XXZ Heisenberg model is introduced in Sec. 4. It 
is used throughout the rest of the paper to fit and interpret 
experimental data. In Sec. 5, the different phases shown in 
Fig. 1 are identified using high-field magnetization data. 
The XXZ model is used to fit the low temperature data 
and to give a physical model for the ordered phase ob- 
served between H& and H c2 . Section 6 is devoted to the 
dynamical processes entering in the NMR relaxation. The 
different regimes decribed in Fig. 1 are presented in Sec. 7 
through T^ 1 measurements across the entire phase dia- 
gram. The Luttinger liquid behavior between H c \ and H c2 
is clearly seen for the first time and compared to the XXZ 
models. Finally the first scaling analysis for a ID quan- 
tum critical point is presented in Sec. 8 and the paper 
concludes with some new perspectives. 



2 Gapped 1-D antiferromagnets: a broad 
universality class 

There are three known families of quasi 1-D materials be- 
longing to the same universality class, with an H-T phase 
diagram similar to Fig. 1. They can all be described by 
the same quantum-field theory, the 0(3) non- linear <j- 
modeljlljnjO] in a magnetic field. Their excitation spectrum |H| 
in a weak magnetic field is represented in Fig. 0. The 
lower critical field H c \ , is reached when the lowest energy 
gap A- vanishes. The upper critical field H c2 is usually 
reached when the highest energy state of the e_ (q) branch 
is below the singlet energy. An approximate representation 
of their singlet ground-states are represented in Fig. [|. 

The know families of 1-D antiferromagnets in this uni- 
versality class are: 

- (i) Quasi lD-antiferromagnetic compounds with two al- 
ternating exchange constants Ji and J 2 have been known 
for over two decades to have a gap (top of Fig. ||). For 
spin- 1/2 alternating chains, the spectrum is identical to 
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Fig. 2. The lowest lying excitations form a triplet branch, 
split by the Zeeman field, or other anisotropic forces. They 
are separated from a unique singlet ground state by energy 
gaps A+, Ao, A-. In absence of anisotropies, they are split by 
the Zeeman energy, A± — Ao ± q^lbH. The shape of the dis- 
persion relation across the Brillouin zone depends on the ex- 
change interactions specific to each 1-D system. In particular, 
the energy gap may be at qa — ir (S=l chains, ladders) or at 
q = (dimerized chains with sufficiently large antiferromag- 
netic second-neighbor exchange). Similarly, a two-"magnons" 
continuum may intersect the triplet branch and have a lower 
energy in some region of the Brillouin zone. In this universality 
class, there are no singlet states in the gap. 



Fig. and in the strong coupling limit (Ji ^> J2), the 

energy gap is A w Jt — 2 ~ §Ji' Very n i ce thermody- 
namic studies[[l5| of Cu(N 03)2.2. 5^0 have identified the 
existence of two critical fields H& ~ Ji — J%/2 ~ 2.8 T 
and H c2 = 2Ji ~ 4.5 T. In spite of the modest value of 
H c i and H C 2, dynamical properties close to these critical 
points have never been thoroughly mapped out neither by 
NMR relaxation measurements nor by neutron scattering. 
Considering the interest in quantum phase-transition, this 
interesting compound should be revisited. 

- (ii) Spin-1 Heisenberg antiferromagnetic chains Q with a 
sufficiently weak planar anisotropy have in zero magnetic 
field a triplet excitation branch separated by an energy 
gap A w 0.41 x J from the unique singlet ground state. 
The most thoroughly studied compound in this family is 
NENP (l6) . Because of the presence of a planar anisotropy, 
this system has three different lower critical fields (9.8, 
13.3 and 14 Tesla) [ p~7| , [T^ ] depending on the orientation of 
the field with respect to the anisotropy axes. The upper 
critical field which has not been measured should exceed 
86 Tesla. Thermodynamic and dynamical measurements Jl9| 
have been carried out at the lower critical field and pro- 
vide very valuable insight on zero-temperature phase tran- 
sitions. 

- (iii) Spin-ladders are quasi-lD structures where a fi- 
nite number of antiferromagnetic chains are coupled by 
a transverse antiferromagnetic exchange. Ladders with an 
odd-number of coupled chains are gapless and belong to 
the same universality class as the spin- 1/2 Heisenberg 
chain pOfl . On the other hand, spin- 1/2 ladders with an 
even number of legs are gapped and form a singlet ground- 
state with short ranged spin correlations (spin- liquid) . While 
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Fig. 3. Three families of gapped antiferromagnets. Top, dimer- 
ized chains: the antiferromagnetic bond Ji dominates J2. Sin- 
glet valence bonds are preferentially localized on these bonds. 
Middle, integer spin chains. The ground state can also be 
viewed as a product of valence bonds [Q, by decomposing each 
spin 1 in two spin- 1/2 and forming valence bonds with the pair 
of spin- 1/2 at the extremity of each bonds (valence bond solid). 
Bottom: spin ladders. In the strong coupling limit (J± 2> J\\), 
singlet bonds are localized on the rungs of the ladder. 



several compounds with ladder- like magnetic structure ex- 
ist |pl| , |22| , the only system where the quantum critical 
point H c i is experimentally accessible is CuHpCl, a coor- 
dination compound made up by stacking binuclear molecules 
in a ladder structure S). Thermodynamic p|, and dy- 
namic quantities ||^|] have been measured over the entire 
phase diagram and give a relatively complete experimental 
picture of a zero-temperature phase transition. 

We now describe its structure, and the relevant mag- 
netic interactions in this material. 

3 CuHpCl, a 1-D spin ladder in the strong 
coupling limit 

The molecular unit is a binuclear structure with two Cu ++ 
ions (spin 1/2) each lying in a middle of two parallel dis- 



torted square-structures) 23 . At the vertices of each square, 
one finds two chlorine and two nitrogen ions as depicted 
in the top of Fig. ||. There are two super-exchange paths 
through the chlorine ions Cl\ and CI2 bridging the cop- 
per ions. However, the p z orbitals of the chlorine ions are 
nearly orthogonal to the Cu ++ d xy orbital. The resulting 
exchange constant between the copper ions J± ~ 13.5 K, 
is found to be weaker than in other materials with similar 
Cu — Cu distance. The organic rings, on the outside of the 
ionic-core just described, contribute further to the delocal- 
ization of the unpaired Cu ++ orbital. Each molecular unit 
stacks up in the [101] direction of this P2\/c molecular 
crystal (see Fig. In addition to the van-der-Vaals forces, 
the organic rings allow a weak hydrogen bonding between 
molecules, along which a super-exchange path can propa- 
gate. In spite of the relatively large intermolecular Cu—Cu 
distance, the orbital overlaps are favorable and lead to an 
intermolecular-exchange interaction Jn w 2A5K along the 
ladder. Strictly speaking, due to the low symmetry, there 



4 Hamiltonian representations of 
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Fig. 4. Top: Representation of the parallel plane geometry of 
the Cu ++ orbitals. Bottom: stacking of the molecular units in 
a ladder structure along the [101] axis. 



are small differences in Cu — Cu distances along the ladder 
above and below each molecular unit. Because the bridg- 
ing entities have the same symmetries, these small differ- 
ences are not expected to modulate significantly the ex- 
change (Jy) along the legs. The translational invariance is 
nevertheless naturally broken, making any additional lat- 
tice (spin-Peierls) instabilities less favorable. Comparison 
of thermodynamic properties to numerical simulations plj 
have demonstrated that the other possible magnetic cross- 
bondings between the legs are weak and need not be con- 
sidered. On the other hand, it is now clear that there is also 
a weak inter-ladder super-exchange, probably also medi- 
ated by a weak hydrogen bonding between organic rings. 
If it is not relevant in the gapped phases, it induces at 
low-T a 3D-ordered phase in the gapless region between 
H c i and ^J C 2 ]3C|] . While this phase has been observed in 
specific heat measurements, its actual structure has not 
been determined experimentally. 



For most purposes, it will be sufficient to consider a quasi- 
ID ladder Hamiltonian in a magnetic field H — H1+H2 + 
Hz, where 

N 

Hi = J± ^2 S 2 i-i • S 2 i (1) 

2N 

H% = J|| ^2 Si ■ S l+ 2, (2) 

1=1 

(even spins are on one leg and odd spins are on the other) . 
The weak ^-factor anisotropy (g± ~ 2.03, <7|| « 2.11) ob- 
served in EPR measurements [p] may be retained in the 
Zeeman Hamiltonian, Hz — J2i=i a 9a^BSfH a . In the 
strong coupling limit (J± ^> Jm), it is possible to give a 
straightforward description of the low-energy states in a 
magnetic field, treating Hi as a perturbation[||[|2^]. The 
eigenstates of H\ + Hz which describes isolated dimers 
in a magnetic field, are built from the singlet J (valence 
bonding) and triplets T +1 , J , (antibonding) on each 
rung. Since we are interested in the critical region, where 
the Zeeman energy is of the order of Jj_, it is legitimate 
to project H2 on the restricted Hilbert space generated by 
the lowest dimers states, J and !_]_. The matrix elements 
of H2 between neighboring dimers can be represented on 
this subspace by a 2 x 2 matrix, which is expressed in 
second-quantized notation as 

J N 

Hf = yJL + *»4f 1 + nrHr+l) • (3) 

r=l 

The fermionic operator v creates the triplet state ?_ 1 on 
bond r = [2i, 2i + l] (only one triplet per bond is allowed), 
while t destroys a triplet, leaving a singlet on bond r. The 
operator n r = t' r t r counts the triplet occupation of bond 
r. In the fermion language, the first two terms represent 
the kinetic energy while the last term is a short range 
repulsion between fermions. Since S 2 and S z are good 
quantum numbers (S is the total spin), it is convenient 
to divide the Hilbert space into sectors with a given value 
of S z . In the restricted Hilbert space, each sector specifies 
the total fermionic occupation since 

S z = J2^r- (4) 

r 

The singlet sector is not coupled by H^ and the singlet 
eigenstate remains the dimer product \S = 0) = |JJ .. .) 
with energy Eq — —jNJ±. This ground state energy 
may be compared to a serie expansion^] in J||/Jj_p7|. For 

1 When the full Hilbert-space is retained, the strong coupling 
expansion for the singlet ground state reads 

\s = o) = \ll-) + ^J2 ^ 
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the parameters appropriate to CuHpCl, the correspond- 
ing singlet energy is only 1.6% lower than in the previous 
estimate. The reduction of H2 to Eq. || is therefore appro- 
priate to CuHpCl, at least for qualitative answers. 

In the S = 1 (1-fermion) sector, the Hamiltonian TL?^ 
can also be diagonalized in Fourier space, since the inter- 
action term n r n r +\ does not contribute. The dispersion 
relation of this triplet state is, 



E Sz= -i{q) = J± - gz^EiH + J|| cos(qa), 

\l,q) =^ eX P(^ r )H---Ir-l*rIr 



(7) 
(8) 



Ao = J i — J11, there is an 



For fields below g z fisH c 

energy gap A- = Aq — g z fJ.sH between the singlet and 
the q = Tt/a minimum of the triplet branch, as shown in 
Fig. g. For fields above H c i or temperatures above the gap 
A-/kB, it is necessary to explore the energy spectrum of 
Ti.2 at finite fermion density. In the low-energy sector of 
the Hilbcrt space (J, a finite fermion density raises 

the energy of Hi + H z with respect to the ground state by 
an amount proportional to the triplet (fermion) density 



6E = (J x - g zf i B H)J2'> 



(9) 



In other words, [i = J± — g z ^BH acts as the chemical 
potential for the triplets. 

The low energy spectrum of Ti^ can be described ex- 
actly on the restricted Hilbert space: Eq. can be rec- 
ognized as the fermion representation of the S = o XXZ 
Heisenberg model 

7Y e ff = J\\ ^ ' ( S r S r _^i + S^S^_^_i + —S r S r _^_i 



(10) 



in an effective field 



H c s = Jj_ + — - g z HBH- 



(11) 



H c ff is zero at the midpoint between H c \ = J± — Ju and 
the upper critical field H c2 = J± + 2J||||. The spin eigen- 
states of this effective Hamiltonian are a representation of 
the triplet-singlet subspace (| f) = *_ l5 1 J.) = J) on each 
rung and have nothing to do with the original spin-1/2. In 
this model, the excitations which carry h/2 angular mo- 
mentum (spinons), have a semionic character, i.e. can only 
be observed in pairs. In non-zero effective fields (H c g 7^ 0) 
the longitudinal and transverse excitations must be dis- 
tinguished. The continuous spectrum of transverse exci- 



( 



+ X 



E = --NJi 
4 



V2 
, 1 ( J H 



.) + 



+ 



(6) 



where the same valence bond notation is used for all singlets 
whether they lie along the legs or the rungs. 




Fig. 5. Energy spectrum of transverse excitations in a mag- 
netic field. Top panel: H > H c i, the incommensurate wavevec- 
tor is close to the antiferromagnetic point. Lower panel: H ~ 
(H cl + H C 2)/2, the incommensurate wave-vector is close the 
zone center. On both plots, the spin-stiffness p a oc sin(7r — Q a ) 
is the slope of the bottom edge of the spectrum at q = n. 
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soliton antisoliton 

Fig. 6. The k = n/2 (a) and k = n (c) ground states are 
superposition of the singlet-triplet configurations shown. The 
wavevector k refers to the wavevector in the XXZ representa- 
tion. In ladder representation shown, this wavevector is fe/2. 
Excitations are built as Bloch wave of the soliton-antisoliton 
pairs represented on (b) and (d). 



tations (1-magnon or 2-spinons) is represented on Fig. |^ 
at two different magnetic fields. At fields just above H cl 
(Fig. U, top panel) a new minimum in the spectrum de- 
velops at Q a a — ira, where a = M{H e g)/M sa t is the 
spin-polarization of the XXZ-model in an effective field 
H c ff. The lower critical field H « H c \ correspond to the 



saturation field of the XXZ model. 



tj sat 



3J|i/2: the 
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spin-polarization a in -ff^ff* is —1 and the ladder magne- 
tization m is zero. Just above H c \, the soft modes at n/a 
and Q a are very close: the spin stiffness p s cx (1 + a) oc m 
goes to zero and a large low-energy spectral weight ex- 
ists close to the antiferromagnetic point. The situation is 
quite different close to H c g — (i.e. half-way between H c i 
and H C 2), where the incommensurate wave- vector Q a , is 
close to the zone center (Fig. [|, lower panel) . For longitu- 
dinal fluctuations, the incommensurate minima in the two 
cases considered are essentially interchanged with respect 
to transverse fluctuations. 

It is useful to "translate" the XXZ incommensurate 
states just described into the valence bond representation 
of ladder states. The pictorial images of the incommensu- 
rate ground states shown in Fig. ^, close to J± — Jn/4 and 
J_l + Jy/2, are appropriate on short lengthscales (quantum 
fluctuations destroy the periodicity Q a on long length- 
scales). Low energy excitation above the ground state are 
Bloch waves of the soliton-antisoliton defects depicted on 
Fig. ^: this builds a coherent superposition aj + /3j_i 
(transverse fluctuation) at wavevector Q a . 

Since so many exact results are known for the XXZ 
model, the mapping discussed here will prove to be ex- 
tremely useful in the analysis experimental data. The high 
field magnetization of CuHpCl clearly establishes the cor- 
respondence with the phase diagram show in Fig. [I]. 

5 Identification of the different phases with 
high-field magnetization measurements 

The weak g-factor anisotropy and the monoclinic sym- 
metry of this crystal, allow a straightforward determina- 
tion of the magnetization by torque magnetometry: if the 
field is applied along the z axis, which does not coincide 
with the principal axes a, b, c, the magnetization is not 
collinear with H, and a torque rocMxfl can be mea- 
sured. The magnetization curves of a 100 fig monocrys- 
tal have been measured with an ultrasensitive AC torque 
magnetometer p8|, p[ . It is straightforward to identify the 
critical field H c \ and H C 2 from the T — 0A2K magnetiza- 
tion curve shown in Fig. ^. Below H c %, the magnetization 
is thermally activated with an effective gap A- (this is 
shown in Fig. |8|). Similarly, above H C 2, the deviation from 
the saturated magnetization are activated with an energy 
gap g z fiB{ii — H C 2) (see NMR section). In appendix A, 
the energy spectrum of excitations carrying one unit less 
angular momentum than the fully polarized state are de- 
termined. In a ladder, there are two spins in the unit cell, 
and hence two spin- wave modes. Since the polarized phase 
is unstable when the lowest energy of these spin-waves 
drops below the energy of the polarized state, the upper 
critical field can be specified exactly 

g z fi B H c2 = J ± + 2J t (12) 

The observed behavior of the magnetization in the differ- 
ent field regions coincides precisely with the zero temper- 
ature phases specified on the y-axis of Fig. 1. Quantita- 
tively, the exchange parameters J± w 13.5 K, Jii « 2.45 AT 




i i i i i i i i i i i i i i i i i i i I 

5 10 15 20 
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Fig. 7. Magnetization of CuHpCl between and 20 T, at 
different temperatures. The symbols are experimental data 
points, while the solid lines are the fits to the XXZ model de- 
scribed in the text, using the procedure outlined in appendix 
B. From the lowest temperature curve, two critical fields H c \ 
and H C 2 are identified. Below H c \, the system is unmagneti- 
zable, i.e. forms a singlet ground state separated from S 7^ 
excitations. Above H c 2, the system is fully polarized. The "in- 
commensurate phase" is between H c i and H C 2- 



are most accurately determined from the values of H c \ 
and H C 2- J± can be identified independently as the mag- 
netic field at which the NMR relaxation rate is maxi- 
mum in the high temperature limit (see Sec. 7). These 
numbers have also been checked against (a) the low and 
high temperature dependence of the susceptibility |^,^9|, 
(b) the gap suppression of the low temperature specific 
heat and (c) the overall bandwidth (2Jii) of the triplet 
branch (H=0)|30| measured by neutron scattering. In a 
numerical study of the ladder magnetization, the pres- 
ence of a weak cross-exchange coupling between legs has 
also been investigated p4| . The conclusion is that this cou- 
pling is weak and if non-zero, ferromagnetic. Considering 
all the experimental and numerical uncertainties, it seems 
at present unnecessary to keep any additional exchange 
coupling in the analysis. 

The thermodynamic properties for the XXZ model can 
be computed exactly by the Bethe Ansatz]3^,^3| and the 
magnetization curves have been evaluated numerically us- 
ing the procedure described in Appendix B. At temper- 
ature below J_l,|35| the results (Fig. [?], solid lines) agree 
very well with the experimental magnetization. Consider- 
ing that there are no adjustable parameters, the mapping 
of strongly- coupled spin-ladders onto the XXZ model ap- 
pears to be excellent. In particular, the quantum critical 
behavior of ladders at H c \ is acurately reproduced by this 
model, a point which is emphasized further in Sec. 8 by 
constructing explicitly the scaling plots for the magneti- 
zation. In the gapped phase, the XXZ mapping progres- 
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H ef f=° H =( H c, + H c2)/ 2 
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Fig. 8. Magnetization of CuHpCl in the gapped phase as a 
function of 1/T, at different magnetic fields. The experimental 
points are compared to the free fermion model of Troyer (solid 
lines). A systematic deviation from the free-fermion model is 
observable when their density becomes significant: this empha- 
size the role of the short range repulsion between fermions. 



sively loses its validity at small Zeeman splitting (com- 
pared to Z\) [^5| . In this limit, it is instructive to compare 
the temperature dependent magnetization to the free-fer- 
mion model proposed by Troyer et al.|34,p| where 



M 



2z(/3) sinh((3g zf i B H) 
1 + [2 cosh (3 g z ^ B H]z{[3) 



(13) 



and z{j3) ss exp(— f3J±)Io((3J\\) In this model, the statisti- 
cal weights are adjusted in order to reproduce the full spin 
entropy at high temperature, but the nearest-neighbor re- 
pulsion is ignored. Fig|| shows that, at high temperature, 
the experimental magnetization is systematically higher 
that inferred by the free-fermion model (solid lines). At 
finite temperature, interactions between fermions increase 
the value of the chemical potential (which is negative in 
the gapped phase, /i > — |Zl_|). Hence at higher tempera- 
ture, the effective gap becomes smaller and a higher 
fermion density is possible. The role of interactions will be 
further emphasized with the identification of the NMR re- 
laxation processes: the dominant relaxation channel (stag- 
gered process) when H — > H cl would not be possible with- 
out interactions. 



5.1 Ordered phase 

Down to T — O.IK, magnetization curves show no plateaux 
nor slope changes which could indicate a 3D ordering tran- 
sition. The ID models appear to give a precise account of 
the magnetization at all temperature. On the other hand, 
small but sharp peaks in the specific heat have been ob- 
served above H c \ at low temperature 1 30, 3l[|. After inte- 
gration of the specific heat at constant field, it is found 
that the entropy per spin associated to this transition is 
very small (< 10 _2 fcB/spin). In light of these experimen- 
tal facts, a second order phase transition appears at low-T 
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Fig. 9. Suggested ordered dimer-state close to half-filling: 
triplet and singlet bonds form a 3D antiferromagnetic lattice. 



between H c \ and H C 2, involving a small change in spin- 
entropy and no detectable change in magnetization. 

When 3D coupling are ignored, the XXZ mapping (Sec. 
4) gives a representation in terms of an interacting ID 
fluid of spinless fermions (Luttinger liquid). Since there is 
no magnetization change, the 3D transition takes place at 
constant fermion density. A very common 3D instability 
for a ID Luttinger liquid is a charge density wave order- 
ing. In a valence-bond language, this transition can be 
viewed as a valence-bond ordering of the ID f i _ 1 states 
in a 3D-lattice. This would hardly affect the magnetiza- 
tion which measures the density of ?_ 1 bonds (fermion 
density) while the quenching of their kinetic energy would 
be manifest in the specific heat. For a 3D-charge order- 
ing, repulsive interactions between fermions are usually 
necessary. Antiferromagnetic superexchange between lad- 
ders, introduce very naturally an additional repulsion be- 
tween fermions on different ladders. This antiferromag- 
netic super-exchange may be represented by, 



Hi 



(k,i) 



(14) 



where the sum (fc, 1} is carried over nearest- neighbor spins 
belonging to different ladders. The physics of a low-T tran- 
sition should be described by the projection of (14) on the 
restricted Hilbert space, i.e. 



rti 



Jin 



E 

(k,i) 



4*i 



to! 



UkUl 



(15) 



which has, up to a sign, the same form as TV^ . A gauge 
transformation, switching the phase of hopping operators 
^2r+i —t2r+i every other rungs, restores an effective 
antiferromagnetic coupling in the spin model. In the spin- 
less fermion model, a transition to a charge density wave 
state can only be established close to half filling [H = 
(H c i + ff e 2)/2l|36|. It is not known whether this tran- 
sition persist at low fermion densities. Since this model 
arise in many context, it will be important to determine 
its complete phase diagram. Of particular relevance is the 
commensurate or incommensurate nature of the 3D charge 
density. The physics of this model is in fact relevant to al- 
most all quasi-lD quantum magnets. For example, in the 
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spin-Peierls compounds CuGeO^, there is also a transi- 
tion to a 3D incomensurate phase, with no change in total 
magnetization. At the transition, only the local distribu- 
tion of magnetization changes [^9t- EBler and Tsvelik|{40f 
have shown that the 3D phonon-couplings present in this 
family of materials can be represented as a transverse ex- 
change between chains in an effective magnetic Hamilto- 
nian. From the point of view of magnetism, this "charge 
density wave" ordering cannot be distinguished from a 
real spin-Peierls transition recently proposed by several 
authors It is therefore natural to expect ordered phases 
with similar structures in all compounds within this uni- 
versality class. 

From the point of view of magnetism, this is an original 
magnetic state, with a 3D ordered structure of valence 
bonds. Fig. |jj gives a pictorial representation of this state 
at half-filling ((H = H c \ + H C 2)/2). From this discussion, 
it is clear that further experimental and theoretical studies 
of the 3D ordering of strongly coupled ladder are called 
for. 



6 Assigment of NMR lines and identification 
of the dynamical relaxation processes in the 
gapped phase 

NMR is an ideal tool to probe the low-energy dynamics 
of quantum magnets ||. When nuclei (here protons) are 
located at different sites than the electronic spins, the in- 
teraction between electronic and nuclear spins are mostly 
dipolar: the dipolar field h.y(i) produced by the electronic 
spin i on the nuclear spin j serves a probe for the dynam- 
ical properties. The time-averaged z-component X^C 1 !?) 
of this local field shifts the value of the magnetic field felt 
by the nucleus by an amount proportional to the local 
electronic susceptibility Xi 



K. 



ijXi 



(16) 



Kj depends on the position of the nuclei in the unit-cell 
through the dipolar sum 



Ai 



_ 2 1 — 3 cos 2 9i j 

-lelnh : ^ - 



(17) 



where Qij is the angle between and H. Kj will be pos- 
itive if hj is mostly parallel to H, and negative if its an- 
tiparallel. CuHpCl contains 24 protons spins in the unit 
cell: in the NMR spectrum at 16 T (i.e. in the polarized 
phase) shown in Fig. [HI there are indeed 24 resolved lines. 
Their position depends on temperature and follow the T- 
dependence of the local magnetization. For 16 nuclei the 
hyperfine shift Kj is positive and negative for the 8 re- 
maining. For proton sites which are nearly equivalents, the 
lines are close together as the local fields (and their fluc- 
tuations) are nearly identical. Using the proton positions 
calculated from X-ray data, it is possible to determine the 
dipolar sum on each site and assign it to a corresponding 
line. For example, the outermost line labeled I arise from 
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Fig. 10. Temperature dependence of the NMR spectrum mea- 
sured at 16T. The intensities have been normalized to compen- 
sate for the gap suppression. Most Tf 1 measurements have 
been carried out on Line I and II, which are always located at 
the extremities of the spectrum for all fields and temperatures. 
They can be identified with nuclei H14&1H23 (line I) and nuclei 
H 2 (line II). 



nuclei H14 & H23 (Ai w 2950 G) while line II arise from 
proton H 2 (An » -2400 G). 

The hyperfine shift is also a measure of the local elec- 
tronic spin susceptibility Mj/H. Its temperature depen- 
dence is shown in Fig. O. Below H c \ there is, at low tem- 
perature, an exponential drop of the hyperfine shift con- 
sistent with an activated behavior with a characteristic 
energy gap Z\_. The dependence observed at 8.7 T, just 
above H c i, shows an increase of the local magnetization 
as the temperature is raised, a very unusual behavior for 
a magnetic phase (e.g. the magnetization of canted anti- 
ferromagnets always decrease with temperature). At low 
temperature (k^T < Jj_) it is meaningful to use the XXZ 
representation, where the system can be viewed as a Lut- 
tinger liquid of triplets: in this limit, the presence of a 
continuum of longitudinal excitation carrying an angular 
momentum h (see Sec. 4) contributes to an increase of the 
triplet occupation with temperature. At higher fields, the 
density of states at small wavevectors gets smaller and a 
more classical behavior is recovered. In quantitative term, 
the hyperfine shift observed below 6K w J±/2 agrees well 
with the XXZ mapping and the thermodynamic measure- 
ments of M/H (see Fig. [y]). 

Temporal fluctuations of the local fields hij (t) at the 
nuclear precession frequency (essentially zero energy) pro- 
vide the dominant relaxation channel for nuclear spins. 
The longitudinal spin- lattice relaxation rate of the nucleus 
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Fig. 11. Temperature dependence of the hyperfine shift at 
different magnetic fields. The thermodynamic measurement of 
Mj H at H=5, 8 and 15T are plotted on the same graph as solid 
lines. The comparison of the NMR shift to the value of M/H 
inferred from the XXZ mapping (Sec. 4) at low temperature 
are drawn as dotted lines. At 8.7 T, in the gapless XXZ-phase, 
there is first an increase of the local magnetization with tem- 
perature, an unusual behavior for a magnetic state. 

is sensitive to the transverse and longitudinal structure 
factors 

S x (q,u n ) = J cxp(-iLu n t)dt(S+(t)Sl q (0)), (20) 
S,(q t u n ) = f exp(-iu n t)dt(S*(t)S z _ q (0)) (21) 

through the form factors F % z and F]_ . These quantities are 
simply the Fourier transform of |^4y | 2 defined in Eq. [l^. 
They are most easily computed in real space as geomet- 
rical dipolar sums: the longitudinal and transverse com- 
ponents of the resultant local field on site i, hj = hy, 
depend on the actual position of the nucleus i in the unit 
cell. Depending on the proton site selected, the relative 
magnitude of the form factors F\ and F]_ can change by 
one order of magnitude: this provides a unique way to 
measure separately all components of the structure factor 
at lu = uj n . 

Fig.[l2| illustrates how different the relaxation on differ- 
ent proton sites can be in the gapped phase (H = 5.6 T). 
If the Tj -1 behavior of both lines I and II are activated, 
the activation energy for line II, An w 6.8 K, is twice as 
large as the activation energy for line I, Aj w 3.4 K. In 
a field of 5.6 T, the smallest activation energy is A- = 
A — g^sH « 3.0 K between the triplet branch |1, q) and 
the ground state, close to the measured energy Aj for line 
I. 

The contribution of one-magnon states |1, q) (cf. Eq. R) 
to the structure factor are proportional to 8(uj — Ei(q)j: 
their spectral weight at the nuclear frequency is zero. On 
the other hand, these states are by no-means exact eigen- 
states and a number of low-energy scattering processes 




Fig. 12. Temperature dependence of the T-[~ relaxation rate 
of sites I and II in the gapped phase (H = 5.6 T). At high 
temperature the relaxation at site II is an order of magnitude 
faster than at site I, but decreases exponentially at lower tem- 
perature with an activation energy An « 6.8 A", twice as large 
as for site I, Ai as 3.4 K. 
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Fig. 13. Schematic representation of multi-magnon scatter- 
ing processes which can drive the NMR relaxation. The in- 
trabranch processes leave m z unchanged and contribute only 
to 5 Z . The interbranch process which involve a spin flip, con- 
tribute to <Sj_ and are quenched one the Zeeman energy (7.6 K 
at 5.6 T) exceed the triplet bandwidth (2Jm ~ 5 K). In the 
staggered process, two magnons from the bottom of the band 

— 1 are scattered into a magnon state at twice the energy: it is 
governed by the square of the magnon density n 2 _ oc exp(2Zi_). 

between magnons can generate a finite spectral weight at 
low energy. Among them, there are (i) finite matrix ele- 
ments of 7^2 between triplet states which are ignored in 
the reduction to the effective Hamiltonian 7^2, (ii) density- 
density interactions (last term in Eq. ||), which contribute 
as the square of the density of thermally excited magnons, 
(iii) other interactions such as interchain couplings and 
impurity scattering. 

It is instructive to follow the classification of low-energy 
processes proposed by Sagi and Affleck j3?j] in the context 
of Haldane S = 1 spin-chains. Only three relevant chan- 
nels need to be examined 

- (i) The simplest processes are the spin-conserving two- 
magnons processes (intrabranch) represented in Fig. [l3| 
the magnons states it — q/2 and 7r + q/2 within a magnon 
branch of a given m z are coupled by the hyperfine interac- 
tion which relaxes the nuclear spins. At low temperature, 
only states in the m z = — 1 branch are thermally occupied 
with a density oc exp(— A^/k^T). This occupation factor 
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dominates the temperature dependence of this relaxation 
channel. In this limit, this process contributes to the lon- 
gitudinal structure factor S(ui n ~ 0, q ~ 0) (no spin-flip). 

- (ii) The spin non-conserving processes (interbranch) cou- 
ple magnons with energies > A in the m z — and m z 

1 branches. They contribute to the transverse structure 
factor (there is a spin- flip), with an activation energy set 
by the full gap A. When the Zeeman energy exceeds the 
one-magnon bandwidth (2J||), which is the case at 5.6 T, 
this process is completely quenched since there are no 
states left in m z = and the m z = — 1 branches with 
the same energies. 

- (iii) Since the magnons states used do not constitute a 
real representation for the eigenstates of the full Hamilto- 
nian, various processes quadratic in the magnon-density 
can occur. The lowest order process contributing to the 
transverse structure factor is a three-magnon process where 
two-occupied magnons at the bottom of the m z = —1 
band are scattered into a magnon with twice the energy 
via a large momentum transfer. The extra angular momen- 
tum being absorbed by the nuclear spin, this process con- 
tributes to the transverse structure factor S±(uj n w 0,q), 
where q is large and will be taken as ir in the rest of the 
analysis. This process is governed by a quadratic thermal 
occupation factor oc n 2 = exp(— 2Z\_//cbT) and requires 
to have a final state in the m z = —1 branch available 
at energy 2Z\_: this is the case when the bandwidth 2Jn 
exceeds the gap A—, i.e. sufficiently close to H c \. There 
are other relevant quadratic processes: four magnons scat- 
tering (Eq. H) processes at the bottom of the m z = — 1 
band have the same temperature dependence but are spin- 
conserving and hence enter only in the longitudinal struc- 
ture factor. 

To summarize, the dominant processes in an interme- 
diate field range are 

- for the longitudinal structure factor S z , the intrabranch 
two-magnon process, S z oc n_ oc exp(— A^/k^T) 

-for the transverse structure factor St, the 3-magnons 
staggered process represented in Fig. S± oc n 2 _ oc 
exp(-2z}_/fc B T). 

In this simple picture, two numbers S z {q = 0) and 
iSj_(<? ~ 7r) are sufficient to specify the T^ 1 relaxation of 
all lines at a given field and temperature 



1 



oc F*{0)S 9 + Fi(ir)Si 



(22) 



Table 1. The form factors for lines I and II, computed as dipo- 
lar sums, are expressed in units of 10 -4 A -6 . The uncertainties 
which are of the order of 20%-30% have a number of origins 
and are discussed elswherefcfl . 
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Fig. 14. Left: activated behavior of T -1 (line I) in the high 
field phase. Right: effective gap as a function of magnetic field. 
The intercept on the field axis coincide with the value of H C 2 = 
13.2 T. 



an order of magnitude. Since the fermion density is like- 
wise large above the critical field H c \ , staggered processes 
dominate the relaxation in the XXZ phase, a result which 
is consistent with all theoretical analysis ]37|, [Tol,!^ and the 
data presented in Fig. [ll| 

The relaxation rate T _1 in the polarized phase (above 
H C 2) is also activated as shown in Fig.pTl The measured 
activation energy is close to g[iB (H — H&) in good agree- 
ment with the exact low energy spectrum described in ap- 
pendix A. Hence, this polarized phase has no Goldstone 
mode and in this sense is not a ferromagnet. There are 
two magnons modes above the polarized state (two spins 
per unit cell): hence, the structure of the spectrum above 
H C 2 and below H c i is qualitatively different, indicating 
that the XXZ model which reproduces correctly the low- 
energy behavior in the gapless phase, cannot be taken lit- 
erally over the entire phase diagram. 



7 Incommensurate phase and quantum 
critical regime 



The value of the form factors F appropriate for line I and 
II are given in Table 1. Since Fj 1 ^ F/ 7 , line II is dom- 
inated by the transverse structure factor, i.e. is governed 
by the staggered processes oc exp(2zi_/fcBT) in this field 
range: this is fully consistent the observed activation en- 
ergy. While longitudinal and transverse fluctuations con- 
tribute to line I, the linear system ( p2| ) can be solved ex- 
plicitly and the different component of the structure fac- 
tor, plotted in Fig. [l^, show very clearly the two different 
activation energies, corresponding to the two dominant 
processes. 

At high temperature, when the thermally excited fermion 
density is important, the staggered processes dominate by 



Fig. 15 gives the overall temperature and field dependence 
of the NMR relaxation rate (line I) through the entire 
phase diagram. In the left panel, the dependence of the 
T^ 1 rate through the lower critical field is displayed. In 
the lower part, the behavior in the gapped phase is re- 
produced over a broader temperature range. Two distinct 
regimes can immediately be recognized. When the tem- 
perature is raised at constant field through the effective 
gap A-, the exponential suppression of the T^ 1 crosses 
over to the temperature independent value of the relax- 
ation rate expected in a classical system. Comparing with 
Fig. |l|, this crossover can be recognized as line (A) (de- 
fined as k^,T = A—) separating the gapped phase, con- 
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Fig. 15. Left panel : temperature dependence of 1/Ti through 
the critical field H c \. In the lower part, the T^ 1 dependence in 
the singlet gapped phase (H < H c \) is displayed, while its be- 
havior in the magnetic phase (H > H c i) is shown in the upper 
part. The dashed lines are guides to the eyes. Lines labeled (A) 
and (C) correspond to the crossover lines of Fig. |l|. Right panel: 
temperature dependence of 1/Ti through the critical field H C 2- 
The upper part shows the behavior in the magnetic phase be- 
low Hci while the lower part is in the "fully polarized" gapped 
phase (H > Hci)- 



trolled by quantum fluctuation, and the quantum critical 
phase where £t is the only relevant lengthscale. When the 
field is raised above H c i, the behavior at high tempera- 
ture is qualitatively the same but at low temperature T^ 1 
starts to diverge. Again, it is possible to place a crossover 
line (C) separating the two regimes. At high temperature, 
one recognizes the same quantum critical phase controlled 
solely by thermal fluctuation, while at low temperature, 
the spectrum of low energy fluctuation in the gapless mag- 
netic phase controls the NMR relaxation. When the mag- 
netic field is raised and crosses H c2 (Fig. |ll| right panel), 
the same features are qualitatively observed, with possi- 
bly a weaker divergence of T^ 1 at low temperature in the 
gapless phase. What is the origin of this critical behavior 
of T^ 1 throughout the gapless phase? 

Because of the broken rotational symmetry, transverse 
(<Sj_) and longitudinal fluctuations (S z ) involve different 
processes (cf. Sec. 6) which do not have the same temper- 
ature dependence. Since staggered processes, entering <Sj_, 
were found to dominate the T^ 1 just below H c i, it is natu- 
ral to first examine the transverse low energy modes in the 
gapless phase. In the XXZ mapping for strongly coupled 
ladders (Sec. 4), two soft modes (transverse in the valence 
bond |, ]|_x basis) were found (Fig. |). One mode is always 



at Q = 7r, and generates the staggered process which was 
already found to be strongly relevant. The other mode is 
at an incommensurate wavevector Q a — ira, which is near 
7T when H is close to H c i (top panel, Fig. ||). In fic|] , it was 
argued that this mode was gapped and did not contribute 
to transverse spin-spin correlation, in apparent contradic- 
tion with the spectrum of the XXZ model discussed in Sec. 
4. On the other hand, this incommensurate soft mode in 
this model is a transverse singlet-triplet wave. For physi- 
cally obscure reasons, the transverse spin-spin correlator 
in this state is indeed found to have zero spectral weight 
at the incommensurate wavevector. This point is crucial 
since it introduces a fundamental difference between spin- 
ladders and integer spin-chains which otherwise belong to 
the same universality class. This has also important con- 
sequences for neutron-scattering studies of spin-ladders. 
For NMR, the staggered process (Q = n) becomes the 
only relevant soft mode for (transverse) relaxation. If the 
temperature exceeds the maximum of the lower edge of the 
spectrum^] between Q a — iro and n (Fig. ^|, top panel), 
many additional modes contribute to S± . Hence the Q = tt 
soft mode is only relevant for temperature below the spin- 
stiffness constant p s , 



7T M 

2 J|I M^ 



J\\9^b(H - H cl ), 



(23) 



where the approximation is valid close to H c \. The con- 
dition kT ps p s specifies the crossover line between the 
'Luttinger liquid' and the quantum critical regime (Fig. [l] 
and [l5| , line C). As the temperature is lowered, the spec- 
tral weight in the Q = tt soft mode increases giving the 
divergent contribution to 



1 



knT 



(24) 



where r\ is the exponent gov erning the power law decay of 
the correlation functions|flC[|. At H c i and H C 2, the expo- 
nent T] is known to be rj = 1/2 [Q and depends smoothly 
on the magnetization in between. The precise dependence 
of r\ with H is not known for ladders and could differ from 
Haldane systems p£|, where the incommensurate mode at 
Q a is relevant. The experimental T-j - divergence observed 
at low temperature appears to be fully consistent with a 
square-root singularity (see Fig. 18). While there could 



also be additional critical fluctuations associated to the 
3D ordering!] occurring at very low temperature, the diver- 
gence of Tj~ which is clearly noticeable below 5K (= 2 J||) 
has to involve ID fluctuations. 

At temperature above p s , the relaxation rate gradually 
crosses over to a constant. This high temperature limit of 
the relaxation rate has a maximum around 10 T, field at 
which g/isH « J±. 

2 This quantity is proportional the spin-stiffness of the anti- 
ferromagnetic magnons. 

3 Considering that the Q = tt magnons in the 3D ordered 
structure proposed in Sec. 4 are essentially the same modes as 
in the ID quantum disordered phase, the ordering should not 
have a dramatic effect on the divergence of the NMR rate. 
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Fig. 16. Field dependence of the longitudinal relaxation rate 
Ti — 1 as a function of magnetic field. The maximum at 10T is 
an independent measure of J±. 



In classical NMR theory jy] , the high temperature limit 
of the T _1 relaxation rate is proportional to the second 
moment of the spectral distribution of excited states. This 
quantity is peaked precisely in the middle of the triplet 
band [m z = — 1). At high temperature where dimers are 
decorrelated, it is thereore natural to find a maximum in 
the relaxation rate as a function of field (shown in Fig. |l6|) 
at the level crossing between the ?_ 1 and J states. Since 
this "classical" contribution is proportional to the zero 
frequency spectral weight, its temperature dependence is 
expected to be weak. It is anyway unrelated to the hydro- 
dynamic soft mode at Q = ir, which is a manifestation of 
the quasi long-range correlation along the ladder. 

The longitudinal fluctuations (<S Z ) have also a contri- 
bution to the relaxation rate originating from the 
Q = uniform mode. They have been shown to be non- 
critical Hoi 



i 



Pa 



(25) 



Since this contribution is noticeable only close to H c \, 
where it is weak (see Fig. |l5|), it will not be discussed 
further. 

Since in the experimental data shown in Fig. [l^, we 
are clearly able to identify the scaling parameters cc< = 
A-/k B T for H < H cl and x> = p s /k B for H > H cl (line 
A and C) on each side of the critical field, it is natural to 
construct the scaling plots appropriate to this quantum 
critical point. 



8 Scaling plots in the quantum critical regime 

The concept of scaling at a quantum phase transition, one 
of the most beautiful idea in condensed matter physics, 
was developed in the context of the metal-insulator tran- 
sition in disordered systems [Q, where it has been bril- 
liantly applied to doped semiconductors p3{ . But it is in 
two-dimensions that the concept has found the most spec- 
tacular applications: disordered superconducting films Jm], 
[45| go directly from a superconducting to an insulating 
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Fig. 17. Magnetization scaling plot. Below H c i all curves 
collapse on a single curve when plotted in terms of the di- 
mensionless scaling variable (A — gpsH) /k&T . As expected, 
the straight line with slope —1 for x > 1 reproduces the 
exponential suppression of the magnetization with an effec- 
tive gap A- at low temperatures. Above H c i, the scaling is 
also excellent (even in the ordered phase) using the parameter 

p 3 _ ^/gMB(-H"--H"cl)J|| 
kBT kBT 



state through a T = quantum phase transition as a func- 
tion of disorder. Josephson-junction arrays have a field 
tune vortex delocalization transition at a critical fraction 
of the flux quantum f c [fl6| . In a two-dimensional electron 
gas, the transitions between quantum Hall plateaux or to a 
Hall-insulating state are also governed by zero temper- 
ature fixed points |)). More recently, the scaling properties 
of a novel metal-insulator transition[Q in silicon MOS- 
FET's have also been thoroughly investigated. In light of 
this, it is surprising to find so few experimental studies [fl9| 
of quantum phase transitions in one dimension. On the 
other hand, many ID systems have Lorentz invariance: 
this confers unique properties to their quantum critical 
points. In particular, their zero temperature critical be- 
havior can be extended to any temperature via confor- 
mal mapping Q]. This enables to give a precise description 
of the finite temperature "quantum critical regime" dis- 
cussed in the introduction, which is most easily revealed 
by a scaling analysis. 

In the last section, the scaling parameters x> = p s /k B 
and x < — A^/k B T above and below H c \ have been iden- 
tified using NMR. They are also appropriate to scale the 
magnetization curves (shown in Fig. ^): the resulting plot 
is shown in Fig. |l7]. Considering the overall quality of 
the scaling, the variable x < and x> are appropriate to 
this quantum critical point. Furthermore, the 0.42K curve 
which crosses into the ordered phase just above H c \ can 
also be included in this plot. 

A similar scaling analysis of the longitudinal relaxation 
rate is also possible. Below if c i, it is straightforward to 
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scale all the experimental T\— 1 curves in terms of the sin- 
gle parameter shown in Fig. [l8] (lower curve). On 
this plot, the crossover line (A) shown on Fig |l| and [l5] 
reduces to the single point x = 1. For x > 1, the exponen- 
tially activated behavior with energy Zi_ is the straight 
line shown with slope —1. Although all curves scale nicely 
for all x, the quantum critical regime is in principle lim- 
ited to fceT < 2Jii. Above H c x, the appropriate scaling 
function is harder to construct because there are noncriti- 
cal (classical) contributions to T-f 1 which need to be sub- 
tracted. For all data shown on Fig. |l5], these contributions 
were determined so that all curves have, after subtraction, 
the same asymptotic limit at high temperature. The field 
dependent constant C{H) substracted was chosen to co- 
incide with the single high-T value of T^ 1 observed in the 
gapped phase (H < H c i, T — > oo) and is close to the high 
temperature limit of T\ — 1 plotted in Fig. [l6| In this way 
the T — > oo limit coincides with the x< = ir> = point 
of the two scaling functions (above and below H c i). After 
subtraction, all data for H > H cl also collapse on a unique 
scaling curve when expressed in terms of the scaling pa- 
rameter x> . Data points with x > 1 are in the Luttinger 
liquid regime and a square-root divergence for T^ 1 (line 
shown) agrees well with the experimental data. The pre- 
cise scaling function can in principle be constructed for all 
x from the zero temperature dynamics]^] and compared 
to this experimental scaling plot. 



9 Conclusions 

In this work, several clean experiments on spin-ladders 
have been used to illustrate the physics of quantum phase 
transitions. Much physical insight could be drawn using 
a mapping of strongly coupled ladders on a much stud- 
ied model of quantum magnetism, the XXZ model. For a 
simple ladder, the anisotropy S = J z /J x = 1/2 puts the 
system in the planar (X-Y) universality class between the 
critical fields H c \ and H c i- Any frustrating coupling (an- 
tiferromagnetic cross-bonding between legs) increases S. 
When 5 > 1, the interactions between fermions (triplets 
bonds are very large and the system switches over to 
an Ising universality class with long range order. When an 
Ising gap is present, plateaux appear in the magnetization 
curve in the vicinity of H c s = (4=> H = (H c i + H C 2)/2, 
half filling) , as the magnetic field has to overcome the Ising 
gap induced by the interactions. This problem has also 
been investigated theoretically by more general means p0|. 
An interesting possibility raised in this paper is that sim- 
ilar physics could be induced by much smaller 3D antifer- 
romagnetic coupling and may already have been observed 
in specific heat experiments [[30|,|3l|] on CuHpCl. 

In the incommensurate gapless phase, low energy dy- 
namical properties are dominated by the Q — n soft mode 
in the transverse excitation spectrum (cf. Fig. |^). These 
fluctuations are not completely quenched in the gapped 
phase, since a three-magnons process, which can easily be 
identified by the temperature dependence of the NMR re- 
laxation rate[^|, remains very strong close to H c \ (T-j -1 oc 



P s /k B T 

,0.0 0.5 1.0 1.5 2.0 
10 i 1 1 1 1 




I i i 

12 3 4 

-A/k B T 

Fig. 18. Scaling plot. Below H c i all the T^ 1 measurements 
collapse on a single curve when plotted in terms of the di- 
mensionless scaling variable (A — q^bH) /ksT. The point x=l 
delimits the quantum critical regime from the quantum dis- 
ordered phase. For x > 1, the straight line with slope — 1 
drawn reproduces precisely the exponential suppression of the 
Tj -1 rate. Above H c i, after subtraction of a constant C(H) 
to the Tj~ rate, it is possible to scale the remaining criti- 
cal contributions to the T^ 1 relaxation in terms of the scal- 
ing variable -jr^Tf = ^— — — 11 . In the Luttinger liquid 
regime (x > 1), the data is consistent with a scaling function 
R(x) — ^Jx (solid line). 



exp(— 2A- /k^T)). At finite temperature, spin-flip pro- 
cesses introduce a natural cutoff of the spin-spin corre- 
lations. When fceT is larger than a characteristic quan- 
tum energy (the effective gap A- below H c \ or the spin- 
stiffness p s above H c i), spin correlation have the same 
nature : this is the quantum critical regime All dy- 
namical properties have in this limit a universal behavior 
which should be common to quantum phase transition in 
ID with a continuous symmetry. The scaling properties of 
the NMR relaxation data on CuHpClQ analyzed in this 
work should serve as a reference for future studies of ID 
quantum critical points. 

In spite of this almost idyllic picture, some questions 
remain open. In the gapless phase, the incommensurate 
soft mode should have observable consequences, probably 
requiring different experimental probes than the ones con- 
sidered here. The nature of the low temperature ordered 
phase appears to hold the answer to several key issues in 
strongly correlated systems. One question to be answered 
is: can ID quantum correlation persist in some form in 
3D ordered phases? Finally, other compounds with more 
frustrated magnetic structures can potentially open new 
horizons in quantum magnetism: for example, there are 
new universality classes which have not been considered 
so far. There are indeed known models|^l| exhibiting an 
energy gap between singlet and triplet sectors but no gaps 
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in the singlet sector. The experimental realization of such 
systems represents a unique challenge in this field. 

The GHMFL is a Laboratoire Conventionne aux Uni- 
versites J Fourier et INPG Grenoble I. 



Appendix A: upper critical field 

The energy of the fully polarized state |/) = | t, 1\ • • •) is 



Er=:,FH F) = x[i± + ^.-g z ^ B H 



The most general state with one unit of angular momen- 
tum less than the completely polarized state is 



1^ = 1(^+^-1)1^ 



(26) 



where the Uj and Vj are respectively the amplitude on the 
lower and upper legs. They are specified by the condition 
that ip) should be an eigenstate of TL with energy , or 
cquivalcntly 



Ep — E- = 



N 

E 



[ujS^+VjS^TiWF). (27) 



This condition yields a set of two coupled equations for 
the Uj and Vj, which are easily solved in Fourier space. 
The dispersion relations for the corresponding spin-wave 
modes e a (k) = E F - E°^(k) are 

e°0) = g z (J, B H- J y (l -cosfc), (28) 
e a (fc) = g z [i B H - J± - Jii(l-cosfc). (29) 



In the strong coupling limit, the 'acoustic' mode (|29|) be- 
comes soft first at wavevector k = tt, when the field drops 
below the field which specifies H C 2 (|l^). It is straightfor- 
ward to generalize the argument to more complex systems. 



the effective field, the thermodynamics is specified by the 
value k(0), i.e. the free-energy per spin is 



Appendix B: thermodynamics of the XXZ model 

This problem has been completely set out by Takahashi 
and Suzuki|33). For an anisotropy factor of J z /J x = 1/2, 
all thermodynamic quantities can be computed from the 
solution of the coupled set of integral equations for the 
functions rj(x) and k(x), 



lnr)(x) = 3V3— + s(x) * ln(l + u(x)) 



k B T 

u(x) = 2k(x) cosh 



2k B T 

\hk(x) = s(x) * ln(l + Tj{x)) 



K 2 {x) 



(30) 
(31) 



where s(x) = ^sech^ and * is the convolution product of 
two functions. For each value of the temperature and of 



£ = I-*bT1hk(0). 



(32) 



We solved Eqs. (p^H) iteratively from the known solutions, 
rj(x) = 3 and k(x) — 2 for J|| = and -ff e // = 0. We 
checked our results against the power series expansion in 

j/tM 



^ = -fc B Tlnf2cosh- 9 ^ 3 
N V 



1 



3J 2 



2k B Tj 8 cosh 2 fjMK. 

V 2k B T 

1/4 



32feT ] cosh 2 juiH.) cosh 4 (^ 



The comparison with the experimental results shown in 
Fig. is really excellent. 
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